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1 Introduction 

The fundamental theory of the strong interaction is quantum chromodynam- 
ics (QCD), which is a relativistic field theory with local gauge invariance, 
whose elementary constituents are colored quarks and giuons. In principle, 
QCD should provide a complete description of nuclear structure and dy- 
namics. Unfortunately, QCD predictions at nuclear length scales with the 
precision of existing (and anticipated) experimental data are not available, 
and this state of affairs will probably persist for some time. Even if it becomes 
possible to use QCD to describe nuclei directly, this description is likely to be 
cumbersome and inefficient, since quarks cluster into hadrons at low energies. 

How can we simplify this problem to make progress? We will employ a 
framework based on Lorentz-covariant, effective quantum field theory and 
density functional theory. Effective field theory (EFT) embodies basic prin- 
ciples that are common to many areas of physics, such as the natural separa- 
tion of length scales in the description of physical phenomena. In EFT, the 
long-range dynamics is included explicitly, while the short-range dynamics is 
parametrized generically; all of the dynamics is constrained by the symme- 
tries of the interaction. When based on a local, Lorentz-invariant Lagrangian 
(density) , EFT is the most general way to parametrize observables consistent 
with the principles of quantum mechanics, special relativity, unitarity, gauge 
invariance, cluster decomposition, microscopic causality, and the required in- 
ternal symmetries. 

Density functional theory (DFT), which has been widely used in atomic 
and condensed-matter physics, allows us to describe the nuclear many-body 
system with a universal energy functional that depends on nuclear densities 
and four-vector currents. In principle, knowledge of the full energy functional 
allows us to calculate any observable for the (zero-temperature) many-body 
system; moreover, a simplified treatment of the functional based on quasi- 
particle orbitals still provides an exact description of bulk nuclear properties 
and some single-particle observables. The great advantage of DFT is that cal- 
culations of this subset of properties can be made without knowledge of the 
many-body wave function, or with a simple one. Finally, if relevant expansion 
parameters can be found, the energy functional can be truncated to a man- 
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ageable size, and the accuracy of the truncation can be tested quantitatively 
for the observables in question. 

The basic properties of nuclei provide stringent constraints on any nuclear 
theory. An accurate description of these properties is necessary for any useful 
predictions or extrapolations. We certainly want to reproduce the observed 
shapes of nuclei: the interior density of a heavy nucleus should be relatively 
constant, there should be a well-defined surface, and because of nuclear "sat- 
uration" , the radius iZ of a nucleus should scale according to i? oc i?^/^, where 
B = N + Z is the total number of neutrons and protons. Moreover, the total 
energy E of the nucleus should agree with the "liquid drop" formula 

E = -aiB + aaB^/^ + a^Z^/B^/^ + ai{N - Zf/B + , (1) 

where typical values for the coefficients are given in ^ |21 El • 

The particle spectrum is determined by the qualitative features of the 
single-particle potential. In a nonrclativistic (Schrodinger) language, the cen- 
tral potential is midway between a harmonic oscillator and a square well; this 
shape determines the ordering of the levels as a function of the orbital angular 
momentum. (See 1 , Figs. 57.1 and 57.2.) In addition, the spin-orbit poten- 
tial is strong, which is instrumental in determining the major shell closures 
and, hence, the nuclear shell model. We will see below how these features are 
easily reproduced in a description based on the Dirac equation. 

These simple nuclear features are the ones we will focus on. We expect 
that they can be adequately described by a single-particle equation with an 
effective, one-body interaction. Such an approach has many names, depending 
on the system being studied and on the practitioner: "shell model" , "mean- 
field theory" , "Kohn-Sham" DFT, etc. Our goal is to correlate (fit) a modest 
number of nuclear bulk and single-particle data and then to predict other, 
similar data as well as possible. 

1.1 Why Use Hadrons? 

Well, why not? Our focus is on low-energy, long-range nuclear characteristics, 
and all measured observables are colorless. (In fact, most of the observables of 
interest to us are dominated by the isoscalar part of the interaction.) More- 
over, hadronic variables (baryons and mesons) are efficient, since hadrons are 
the particles that are observed in experiments. Colored quarks and gluons 
participate only in intermediate states, and such "off-shell behavior" is unob- 
servable; by using hadrons, we expend no theoretical effort combining quarks 
and gluons into color singlets that can actually be observed. 

So we pick the most efficient degrees of freedom by choosing hadrons. We 
will have to parametrize the nuclear EFT Lagrangian anyway, since we cannot 
compute its true form from QCD, and hadronic variables, if combined in all 
forms consistent with the underlying symmetries, provide sufficient flexibility 
for our parametrization. We cannot guarantee that a single-particle hadronic 
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approach will be successful in describing all of the observables of interest, but 
we want to see how well we can do. 

1.2 Why Use the Dirac Equation? 

To motivate the Dirac equation as straightforwardly as possible, compare 
the particle spectrum (and fine structure) in a light atom with the spectrum 
in a heavy nucleus. An example of the former is given in while an early 
example of the latter is given in [H], which is reproduced in Fig. 57.3 of Q]. The 
most striking result is that it is impossible to draw the atomic fine structure 
to scale, since the splittings are roughly 1/10,000 as large as the major- 
level splittings (at least for the deeply bound atomic levels). In contrast, the 
nuclear spectrum shows that the "fine" structure is really "gross" ; the spin- 
orbit splittings are as large as the major-level splittings to within a factor of 
two! 

The implication is that there must be some relativistic effects that are 
important in nuclei (unlike light atoms), and thus it is much more natural to 
use the Dirac equation to describe the quasi-particle nucleon wave functions. 

1.3 Quantum Hadrodynamics (QHD) 

We will refer to Lorentz-covariant, meson-baryon, effective field theories of 
the nuclear many-body problem as "quantum hadrodynamics" or QHD |31 
EllZIIHlElIinillllinilSl- When QHD is applied within the framework of 
modern EFT and DFT, it provides a quantitative description of bulk nuclear 
properties and the spin-orbit force throughout the Periodic Table |H1 
[TT] . This success arises from the presence of large Lorentz scalar and vector 
mean fields, which imply that there are large relativistic interaction effects 
in nuclei under normal conditions There is evidence from QCD sum 
rules that these large fields are dynamical consequences of the underlying 
chromodynamics |15[I16| . Moreover, similar relativistic effects are responsible 
for the efficient description of spin observables in medium-energy proton- 
nucleus scattering using the Relativistic Impulse Approximation, and they 
are consistent with the major role played by scalar and vector meson exchange 
in modern boson-exchange models of the nucleon-nucleon (NN) interaction. 
All of these features motivate further investigation into the application of 
QHD to the nuclear many-body problem. 

2 Effective Field Theory 

A modern discussion of QHD begins by interpreting the Lagrangian as defin- 
ing a nonrenormalizable, effective field theory (EFT) ^| |^. An effective 
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Lagrangian consists of known long-range interactions constrained by sym- 
metries and a complete, non-redundant set of generic short-range interac- 
tions (i.e., "contact" and "gradient" terms). The division between "long" 
and "short" is characterized by the breakdown scale A of the EFT. While it 
is not possible at present to derive an effective hadronic theory directly from 
the underlying QCD, the EFT perspective implies that this is not necessary. 
If one constructs a general Lagrangian that respects the symmetries of QCD: 
Lorentz covariance, parity conservation, time-reversal and charge-conjugation 
invariance, (approximate) isospin symmetry, and spontaneously broken chiral 
symmetry, then the EFT is a general paramctrization of observables below 
the breakdown scale. ^ 

For QHD, we identify A with the mass scale of physics beyond the Gold- 
stone bosons (pions); we will see that A « 600 MeV. At momenta small com- 
pared to A, short-distance physics (such as the substructure of nucleons) is 
only partially resolved and so may be incorporated into the coefficients of field 
operators organized as a derivative expansion. The coefficients of these short- 
range terms may eventually be derived from QCD, but at present, they must 
be fitted by matching calculated and experimental observables. In principle, 
there are an infinite number of (nonrenormalizable) terms, but in practice, 
the Lagrangian or energy functional can be truncated to work to a given 
precision |H] . The EFT is useful if this truncation can be made at low enough 
order that the number of free parameters is not prohibitive. 

The EFT perspective, with the freedom to redefine and transform fields, 
implies that there are infinitely many representations of low- energy QCD 
physics. But they are not all equally efhcient or physically transparent. One 
of the possible choices is between Lorentz-covariant and nonrelativistic for- 
mulations. Recent developments in baryon chiral perturbation theory support 
the consistency (and utility) of a covariant EFT, with Dirac nucleon fields 
in a Lorentz-invariant, effective Lagrangian density ^1^1120]. A similar 
framework underlies QHD approaches to nuclei. 

In QHD, the only essential degrees of freedom are the nucleons and pions. 
Only these stable particles can appear on external lines with timelike four- 
momenta. The long-range pion-pion and pion-nucleon interactions are in- 
cluded in a nonlinear realization of the spontaneously broken SU(2)^ xSU(2)^ 
chiral symmetry, which avoids dynamical assumptions inherent in linear rep- 
resentations. These interactions can be written down systematically, given 
an appropriate power-counting scheme, to be discussed shortly [S]. Low-mass 
vector mesons are typically included for phenomenological reasons, but are 
not required, since their masses are roughly equal to the breakdown scale 
A; they are absent from point-coupling Lagrangians, for example |^ \22\ . 
In descriptions of NN scattering and of nuclear structure and reactions, the 
heavy, non-Goldstone bosons appear only on internal lines (with spacelike 

^ It is straightforward to include the local U(l) gauge symmetry of the electro- 
magnetic interaction ^_9,. 
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four- moment a) and allow us to parametrize the medium- and short-range 
parts of the NN interaction, as well as the electromagnetic form factors of 
the hadrons |^ |H| ■ The heavy bosons are also convenient degrees of freedom 
for describing nonvanishing expectation values of bilinear nucleon operators, 
such as A^A^ and Nj'^N, which are important in nuclear many-body systems 
Plin]. This explains why it is useful to introduce collective degrees of freedom 
with other quantum numbers, such as a A baryon to incorporate important 
pion-nucleon interactions ^1 ES]- Because one must always truncate the 
EFT Lagrangian, these degrees of freedom can be efficient in the many-body 
problem whether or not they are actually observed as hadronic resonances. 

A Lorentz scalar, isoscalar mean field in nuclei is an efficient way to in- 
clude implicitly the effects of pion exchange that are the most important for 
describing bulk nuclear properties. Because the chiral symmetry is realized 
nonlinearly, one can add to the theory a light scalar, isoscalar, chiral-singlet 
field with a Yukawa coupling to the nucleon, just as in the original Walecka 
model 0. Nonlinear self- interactions of this new scalar must be included, 
with adjustable couplings that arise in part from the nucleon substructure. 
Since the expectation value of the pion field in nuclear matter vanishes at the 
mean-field level, one makes the remarkable observation that the mean-field 
theory (MFT) of the Walecka model is consistent with chiral symmetry, pro- 
vided we think in terms of a nonlinear realization of the symmetry. The light 
scalar, isoscalar field, which is not the chiral partner of the pion, plays the 
same role in the EFT as it does in the Walecka model: it simulates important 
TTTT and NN interactions that must be included from the outset to generate a 
realistic description of nuclear matter and nuclei. 

To make systematic calculations, the EFT approach exploits the separa- 
tion of scales in physical systems, with the ratios of scales providing expansion 
parameters. A connection between appropriate QCD scales and nuclear phe- 
nomenology is made by applying Georgi and Manohar's naive dimensional 
analysis (NDA) ESI and naturalness, namely, that all appropriately de- 
fined, dimensionless couplings are of order unity. With this input, the non- 
linear chiral Lagrangian can be organized in increasing powers of the fields 
and their derivatives. To each interaction term we assign an index 

TL 

y = d + - + b, (2) 

where d is the number of derivatives, n is the number of nucleon fields, and h is 
the number of non-Goldstone boson fields in the interaction term. Derivatives 
on the nucleon fields are not counted in d because they will typically intro- 
duce powers of the nucleon mass M, which will not lead to small expansion 
parameters. 

It was shown in j26[ |H] that for finite-density applications at and below 
nuclear matter equilibrium density, one can truncate the effective Lagrangian 
to terms with < 4. It was also argued that by making suitable definitions 
of the nucleon and meson fields, it is possible to write the Lagrangian in 
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a "canonical" form containing familiar noninteracting terms for all fields, 
Yukawa couplings between the nucleon and meson fields, and nonlinear meson 
interactions |27| ■ See |H1 El for a more complete discussion. 

If we keep terms with < 4, the chirally invariant Lagrangian can be 
written as^'^ 

-Ceft = -Cat + £4 + 

= N{i-f'' [d^ + iv^ + igpp^ + ig^Vf,] + 5^7^750/. - M + gs(f>) N 

- '^Nvp^a^^'^N + ^NN Tr {a^a^ + £4 



M ^ 

+ ^ Q^gl{V,Vn^ + (1 + ml Tr {p,p^) 

where the nucleon, pion, sigma, omega, and rho fields are denoted by N , tt, 4>, 
Vf,, and p^ = ^T-p^^, respectively, Vf,^ = d^,V^ - d^V^, and cr'"' = 
The trace "Tr" is in the 2x2 isospin space. The pion field enters through 
the combinations 

U = exp(iT-7r//^) , ^ = exp(zT-7r /2/^) , (4) 

a^^~\{ed^i-£,d^e)^al, (5) 

y^^^l{S^'^d,(,+^d,e)=vl. (6) 

Vf_,„ = di_,Vu - dyV^ + i[v^, Vy] = -i[a^, a^,] . (7) 

The rho meson enters through the chirally covariant field tensor 

^ We use the conventions of [S1IH1I3- The pion-decay constant is ~ 93MeV. 
^ The two terms involving a.i coefficients actually have = 5, but they were found 
to be numerically significant in See Fig. 1^1 below. 
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Pt^v = D^^pi, - D^^Pf, + i gp[p^, py] , (8) 
where the covariant derivative is defined by 

D^p^ EE df_,p„ + i[vi_,,pu\ ■ (9) 



The antisymmetric combination of derivatives in p^^^ impUes that the 
timehke components of the rho field have no conjugate momenta and 
are thus determined by equations of constraint, as appropriate for a massive 
vector field with three dynamical degrees of freedom. The final term in JHJ 
has the usual form for a nonabelian vector field and enables the rho meson 
to couple to a conserved isovector current |31 121]. C4, contains tttt and ttN 
interactions of order v = A that are not needed in this work. A numerically 
insignificant v — A term proportional to 0^ Tr {ppp'^) has been omitted. 

To exhibit the chiral invariance of £eft explicitly, we follow CCWZ |29l 
131)1 A nonlinear representation of the chiral group SU(2)^ x SU(2)^ is 
defined such that for arbitrary global matrices L G SU(2)^ and R G SU(2)^, 
there is a mapping 

L®R: {^,p^,N)-^{e,p'^,N') . (10) 

Because of the parity operation V, which produces the transformation 

V: L^R, ^-7r«(t,-x) , ^(i, x) ^ e^(^, -x) , (11) 

the chiral mapping (|l()|l can be written as [21] 

('{x) = L£,{x)h\x) = h{x)i{x)R'^ , (12) 

p'pix) = hix)p^,ix)h\x) , (13) 

N'{x) = h{x)N{x) . (14) 

The second equality in (|12() defines h{x) as a function of L, R, and the local 
pion fields: h{x) = h{L, R,tv{x)). It follows from (|12|) that h{x) is invariant 
under the parity operation l|ll|l . that is, 

h{x) e SU(2)^ , (15) 

where SU(2) Y is the unbroken vector subgroup of SU(2) ^ x SU(2)^. Note that 
the matrix h{x) becomes a constant only when L = R, so that h = L = R. 
Equations H13() and H14() then ensure that the rho and nucleon fields transform 
linearly under global SU(2)y, in accordance with their isospins. The isoscalar 
fields V^{x) and (j){x) are chiral scalars and are unaffected by both chiral and 
isospin transformations. 

For discussing purely pionic interactions, it is convenient to use the matrix 
U{x) of Q), since the transformation law (|12|1 then implies 



U{x) — > U'{x) = LU{x)R'' , 



(16) 
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SO that U (x) always transforms globally. Thus derivatives of U (x) transform 
the same way as U{x), and chiraUy invariant interactions involving pions 
afone can be constructed from products of U{x), U^x), and their derivatives. 
As is weU known, these terms can be organized according to the number of 
derivatives, resulting in the Lagrangian of chiral perturbation theory |2 1^ ■ 
We will return to this later when we discuss electroweak interactions with 
nuclei. 

For describing the interactions of pions with other particles, U{x) is not 
convenient, because other fields transform with the local function h(x) of 
the unbroken isovector subgroup SU(2)y. It follows from the transformation 
laws given earlier that interaction terms that are invariant under local isospin 
rotations will be invariant under global transformations of the full group 
SU(2)^ X SU(2)j^. Thus, to form chirally invariant interactions involving pions 
and other fields, we need functions of the pion field that transform with h{x) 
only. 

The desired functions involving one derivative of the pion field are given 
in © and 0. The parity transformation implies that is an axial 
vector and is a polar vector. Moreover, under a chiral transformation, 
(221) implies 

— > a'^ = ha^h) , (17) 
— > v'^ = hv^ih^ — ihdf^h^ = hv^h'' + i{d^ih)h^ . (18) 

Thus transforms homogeneously under the local SU(2)y group and can 
be interpreted as a covariant derivative of the pion- field matrix ^(x). In con- 
trast, the inhomogeneous transformation law for resembles that of a gauge 
field, so that allows us to construct chirally covariant derivatives of the 
other fields. For example, it is straightforward to verify that the covariant 
derivatives @ and 

D^,N = {d^ + iv^,)N (19) 
transform homogeneously with h{x) under the full group: 

{D^N)' = h{D^N) , {D^p,)' = h{D^p,)h^ . (20) 

The covariant tensor for the pion field is [see Q], which transforms 
homogeneously with h, as does p^^. This allows us to produce a chirally 
invariant pmr coupling through an interaction of the form Tr (pi^^^v^^). 

Electromagnetic interactions can be included by adding a chirally nonin- 
variant Lagrangian [we exhibit terms of 0(e) only] 

Ckm = -~\ F^.F'^'' - 2ef^A^ Tr {v^t^) - ^ F^, [Tr (p^Vg) -I- ^V^-"] 



2 

e 



T3 a-ui^N 



2A/P 



d^Ff"" N 



(/3(°)-H/?Wt3)7m]^, (21) 
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where A'^ is the photon field and = d^A^ — d^A^ is the usual Maxwell 
tensor. The constants e, A*^*^ , and in H21|l (where t = 0,1 denotes the 
isospin), together with the tensor couplings /v and fp in are sufficient 
to parametrize the empirical nucleon charge e, the anomalous moments Ap^n, 
and the charge and magnetic radii {rrmB)i \ momentum transfer 

The expansion can be extended to include higher derivatives of the photon 
field if greater accuracy is needed . 

Similarly, the free parameter gA in the pion-nuclcon interaction [see Q] 
allows us to normalize the one-body, axial-vector nuclear current so that the 
Goldberger-Treiman relation is satisfied at the tree level p2j . 

To summarize the important points of the full Lagrangian >Cqhd = 
£eft + -Cem [recall @ and ^] : 

• The noninteracting hadron terms take their standard canonical forms. 

• The generalized coordinates (fields) have been chosen so that the meson- 
nucleon couplings have a simple Yukawa form. 

• The pion-nucleon and pion-meson interactions enforce the nonlinear re- 
alization of chiral symmetry. 

• The nonlinearities involving chiral singlet fields are obviously invariant, 
and fitting their coefficients to data will implicitly include short-range 
dynamics from many-nucleon forces, fluctuations of the quantum vacuum, 
and hadron substructure. 

• The nucleon electromagnetic (and weak) structure ((7A,A,etc.) are in- 
cluded to the desired accuracy using a derivative expansion of the fields. 

3 Density Functional Theory 

The successes of QHD mean-field phenomenology are, at first, rather mys- 
terious from the EFT perspective alone, since the Hartree approximation is 
just the finite-density counterpart of the Born approximation at zero density. 
The density functional theory (DFT) perspective explains the successes of 
mean-field approaches and provides a new context for EFT power counting. 

We begin with a discussion of nonrelativistic DFT and generalize later to 
include relativity. The basic idea behind DFT is to compute the energy E 
of the many-fermion system (or, at finite temperature, the grand potential 
^2) as a functional of the particle density |22|- DFT is therefore a successor 
to Thomas-Fermi theory |84l 1851 . which uses a crude energy functional, but 
eliminates the need to calculate the many-fermion wave function. 

The strategy behind DFT can be seen most easily by working in analogy 
to thermodynamics |36| . For a uniform system in a box of volume V at 
temperature T, one first computes the grand potential f2{fi,T,V), where fj, 
is the chemical potential. It then follows that the number of particles N is 
determined by ^ 
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According to Gibbs, thermodynamic equilibrium is defined by the condition 

{5n)t..T.v > ; (23) 

an assembly minimizes its thermodynamic potential at fixed /i, V, and T. 
Thus the convexity of f2 implies that iV is a monotonically increasing function 
of fi, so the relation H22() can be inverted for /i(iV). Finally, one makes a 
Legendre transformation to the Helmholtz free energy 

F(N,T,V) = n{fi{N),T,V)+pL{N)N (24) 

to discuss systems with a fixed density n = N/V. 

For a self-bound, finite system, we replace the chemical potential with 
an external, single-particle potential^ X^i The grand potential is now 

a functional: J7[w(r);T), and a functional derivative with respect to v gives 
the particle density:'^ 

= {n{r)) - ^ . (25) 

The convexity of f2 allows us to invert this relation (in principle) and to find 
v{r) as a (complicated) functional of n(r): 

v{r) = v[n{r)] . (26) 

(T is suppressed.) Thus there is a one-to-one relation between the external 
potential and the particle density. Needless to say, the possibility of this 
complicated inversion is a matter of great technical interest, and the reader 
is referred to PHI for a discussion. 

Finally, we make a functional Legendre transformation to define the 
Hohenberg-Kohn free energy, which is a functional of n(r): 

FnKHr)] = n [v[n{r)] ] - ^ dr n{r)v{r) . (27) 

The variational derivative of this free-energy functional with respect to n now 
gives 

^ = -(r) , (28) 
dn(r) 

where we have used (|^ . 

If we now restrict consideration to T = and v{r) = 0, the free energy 
becomes simply the energy: -F'HKi?^] — ^ ^'hkM, and the Hohenberg-Kohn 

* Here the chemical potential fi is absorbed into the definition of v, which defines 
the zero of energy. We suppress all spin dependence at this point. 

* Higher variational derivatives yield various correlation functions. 
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theorem follows |Slin7|: If the functional form of i^HKi^-lr)] is known exactly, 
the ground-state expectation value of any observable is a unique functional 
of the exact ground-state density. Moreover, it follows immediately from H28|l 
that the exact ground-state density can be found by minimizing the energy 
functional. Although we have assumed here that the ground state is non- 
degenerate, this assumption can be easily relaxed 

Significant progress in solving these equations was made by Kohn and 
Sham |SH1, who introduced a complete set of single-particle wave functions. 
The exact Hohenberg-Kohn free energy for an inhomogeneous (finite) many- 
body system in an external potential takes the form 

i^HK [n(r)] = ^^„i [n(r)] + E,^, [n(r)] , (29) 

where the subscripts "ni" and "int" denote noninteracting and interacting, 
respectively. -F'ni['^(r)] represents the kinetic energy contribution. The interac- 
tion energy Eii-it[n{vy\ is some functional of the density (and its derivatives); in 
the many-body problem, it contains a Hartree term, an exchange-correlation 
("xc") contribution, etc. |T]: 

Eint[n] = i^Hartrcoi^l] -|- E'xc ["-] H ■ (30) 

Note that E^c is generally a nonlocal and nonanalytic functional of the density 
that contains both many-body and short-distance physics, including vacuum 
fluctuations and hadron substructure. 

Now consider the (nonrelativistic) Schrodinger equation in a potential 
i;off(r), which is designed to give the exact density n(r): 



"2m 



V^-f i;eff(r))^,(r) = e,V^(r) , (31) 

N 

n(r) = 5]|V.(r)r- (32) 

i = 1 

In this problem, VeS plays exactly the same role as the previous u, and thus 
the Hohenberg-Kohn equation H28I) gives 

SFni[n] 

= -eff(r) , (33) 

By taking the variational derivative of (|29|) with respect to n{r), using ()28|l . 
and rearranging terms, we find 

vM . .(r) + ^ . (34) 
on[r) 



Upon setting the external v{r) = 0, we obtain the effective potential to be 
used in the Kohn-Sham (KS) equations H31|l : 
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5n{r) 



(35) 



Thus, if we know the exact interacting energy as a functional of the den- 
sity, we can reproduce the exact interacting density using a set of single- 
particle wave functions. Kohn calls these the "density-optimal" single-particle 
wave functions jSZI (as opposed to Hartree-Fock wave functions, which are 
"total-energy optimal" ) . 

The generalization of DFT to relativistic systems is straightforward [SHI ■ 
The energy E'hk now becomes a functional of both the ground-state scalar 
density ps and the baryon four-current density B^. Extremization of the 
functional gives rise to variational equations that determine the ground-state 
densities Ps and p^ = -Bo- 

These equations can again be simplified by following the Kohn-Sham 
approach. In the relativistic case, the complete set of single-particle wave 
functions allows us to recast the variational equations as Dirac equations 
for occupied orbitals. The single-particle Hamiltonian contains local, density- 
dependent, Lorentz scalar and vector potentials, even when the exact energy 
functional is used. Moreover, one can introduce auxiliary (scalar and vec- 
tor) fields ^{x) and W{x), which correspond to the local potentials and can 
therefore be identified as relativistic KS potentials. The auxiliary fields are 
determined by extremizing the energy functional, which gives rise to a Dirac 
single-particle Hamiltonian. The isoscalar part (for spherical nuclei) looks like 



where M is the nucleon mass and we define M* = M — The resulting 
coupled differential equations resemble those in a relativistic MFT calculation 
IHlini. Note that <I> need not be simply proportional to the isoscalar, scalar 
field (j). In fact, $ could be proportional to (j> (as in the Walecka model), or 
could be expressed as a sum of scalar and vector densities (as in relativistic 
point-coupling theories) , or could be a nonlinear function of i/i (as in modern 
chiral EFT's). 

The strength of the KS approach rests on the following theorem: 

The exact ground-state scalar and vector densities, energy, and chem- 
ical potential for the fully interacting many-fermion system can be 
reproduced by a collection of (quasi) fermions moving in appropriately 
defined, self-consistent, local, classical fields. 

The proof is again straightforward |37| . Start with a collection of noninteract- 
ing fermions moving in an externally specified, local, one-body potential. The 
exact ground state for this system is known: just calculate the lowest-energy 
orbitals and fill them up.^ Therefore, if one can find a suitable local, one- 
body potential based on an exact energy functional, the exact ground state 

^ For simplicity, we assume that the least-bound orbital is completely filled, so the 
ground state is non-degenerate. 



iV ■a + I3[M - $(r)] -f W{r) , 



(36) 
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of that system can be determined. But this potential is precisely analogous 
to i'eff(r) discussed above, which is obtained by differentiating the interaction 
parts of -Fhk with respect to the various densities. The resulting one-body 
potential will generally be density dependent and thus must be determined 
self-consistently. 

Several points are noteworthy. Since the single-particle basis constructed 
as described above is again "density optimal" , the exact scalar and vector 
densities are given by sums over the squares of the Dirac wave functions, with 
unit occupation probability. Moreover, since these densities are guaranteed 
to make the energy functional stationary [the external v{r) = 0], the exact 
ground-state energy is also obtained. The proof that the eigenvalue of the 
least-bound state is exactly the Fermi energy is given in 001 • Note, however, 
that aside from this association, the exact Kohn-Sham wave functions (and 
remaining eigenvalues) have no known, directly observable meaning. 

If one knows the exact functional form of the energy on the densities, 
one can describe the observables noted in the theorem exactly (and easily) 
in terms of the Kohn-Sham basis. Observables of this type are typically the 
ones calculated in relativistic MFT's ^| 021 El- Moreover, it has been known 
for many years |48l El that the mean-field contributions dominate the single- 
particle potentials at ordinary densities. Thus, by parametrizing the energy 
functional in a mean- field (or "factorized" ) form, and by fitting the parame- 
ters to empirical bulk and single-particle nuclear data, one should obtain an 
excellent approximation to the exact energy functional in the relevant density 
regime. This is the key to the success of relativistic MET calculations, as we 
will verify below, using the effective chiral Lagrangian constructed in Sect.O 

4 Naive Dimensional Analysis 

There is still an important point to be addressed: we must understand how 
to extract the dimensional scales of each term in the Lagrangian, so that the 
remaining dimensionless constants can be checked for naturalness. A naive 
dimensional analysis (NDA) for assigning a coefficient of the appropriate size 
to any term in the effective Lagrangian has been proposed by Georgi and 
Manohar |24l 1^ . This allows for a determination of both the dimensional 
scales associated with each term and for the inclusion of an overall dimen- 
sionless constant that can be used to adjust the strength. The basic idea of 
naturalness is that once the appropriate dimensional scales have been ex- 
tracted, the overall dimensionless coefficients should all be of order unity. 
The NDA rules for a given term in the Lagrangian density are [2^: 

1. ) Include a factor of l//,r for each strongly interacting field. 

2. ) Assign an overall factor of /^A^. 

3. ) Multiply by factors of 1/A to achieve dimension (mass)^. 

4. ) Include appropriate counting factors (such as 1/nl for 0"). 
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Here /^r ~ 93MeV is the pion-decay constant, and the breakdown scale 
A « 600 MeV is taken as the generic large-momentum-cutoff scale, which 
characterizes the mass scale of physics beyond the Goldstone bosons. 

As noted by Georgi [53], rule (1) simply assumes that the amplitude for 
producing any strongly interacting particle is proportional to the amplitude 
/tt for emitting a Goldstone boson. This is a reasonable assumption, since /tt 
is the only natural scale. Thus, by dividing each field by /tt, we should arrive 
at a factor of 0(1). Rule (2) can be understood as an overall normalization 
factor that arises from the standard way of writing the mass terms of non- 
Goldstone bosons. For example, one may write the mass term of a isoscalar, 
scalar field Mx) as 



-f'^'Tjll- 

where the scalar mass mg is treated as roughly the same size as A. By applying 
rule (1) and extracting the overall factor of /^A^, the remaining ratios are 
of 0(1). Since all terms will have the same overall scale factor /^A^, higher- 
order terms or terms with gradients of fields will be suppressed by powers of 
1/A relative to the leading mass terms, as a result of "integrating out" physics 
above the scale A. (A simple example is the low-momentum expansion of a 
tree- level propagator for a heavy meson of mass mn, which leads to terms 
with powers of d^/mjj.) It is precisely because of these 1/A suppression 
factors and dimensional analysis that one arrives at rule (3). The origin of 
the combinatorial factors in rule (4) is discussed in 

Applying these rules to a generic term in an eff'ective Lagrangian involving 
the isoscalar fields and the nucleon field leads to (the generalization to include 
the pion, rho, and photon is straightforward) [^15] 

where ip is a baryon field, F is any Dirac matrix, derivatives are denoted 
generically by d, and we have allowed for the possibility of chiral-symmetry- 
violating terms that contain the small parameter mjr/A. The product of all 
the dimensional factors then sets the scale in terms of the pion-decay con- 
stant /tt and the EFT breakdown scale A. The overall coupling constant C 
is dimensionless and of 0(1) if naturalness holds. 

These scaling rules imply that a general potential for the scalar meson 
can be expanded as 

Vs = ^A (^2+3! + ¥ M^ + -"j ' 

in agreement with the corresponding term in Here we have included a 
factor of I/Ztt for each power of 0; these factors are then eliminated in favor of 
gs ~ M / fi^, which is basically the Goldberger-Treiman relation i2j. Factorial 
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counting factors are also included, since the NDA rules are actually meant to 
apply to the tree-level scattering amplitude generated by the corresponding 
vertex 

The naturalness hypothesis states that after the dimensional factors and 
appropriate combinatorial factors are extracted, the overall dimensionless 
coefficients [C in I^Hjl and K3, K4, ... in Ip!?^ ] should be of order unity. It 
should be clear, however, that the preceding arguments are not a proof of 
naturalness, since unknown physical scales could generate unnaturally large 
coefficients. Moreover, some fitted constants may be unnaturally small, which 
often signals a symmetry of the theory that has not yet been identified. These 
caveats notwithstanding, without the naturalness hypothesis it is basically 
impossible to construct an effective Lagrangian with any predictive power. ^ 
Until one can derive the effective hadronic Lagrangian from QCD, the nat- 
uralness hypothesis must be checked by fitting to experimental data, as we 
will do in the following section. 

If truncations of the EFT Lagrangian determined by NDA and naturalness 
are valid, it should also be possible to determine the level of truncation that 
exhausts the information content of the input data. In other words, we should 
be able to verify that adding terms with higher values of v does not improve 
the fits to the empirical observables of interest. To this end, it is useful to 
consider the v — b interactions 

to check that their contributions are either negligible or can be absorbed into 
slight modifications of the parameters in the v < A Lagrangian £qhd . 

5 Mean-Field Theory of Nuclear Structure 

The mean-field equations and energy resulting from /Cqhd in ^ and l|21|l can 
be derived straightforwardly. The interested reader is referred to [H|, where 
the procedures and results are discussed thoroughly. One important result 
is that due to the additional nonrenormalizable interactions in £em between 
the nucleon and the electromagnetic field, and also due to vector-meson domi- 
nance, the computed nuclear charge density automatically contains the effects 
of nucleon structure, and it is unnecessary to introduce an ad hoc form factor. 

As we will show in Sect. El the full MFT Lagrangian derived from >Cqhd 
has more than enough parameters to accurately describe the desired nuclear 
properties discussed in Sect.^ The more important question is whether the 

The assumption of renormalizability also leads to a finite number of parameters 
and well-defined predictions, but does so by imposing unnatural restrictions on 
the Lagrangian, namely, that many parameters are identically zero in the absence 
of relevant symmetry arguments. 
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parameters fitted to nuclei are natural. In 8 , the parameters were determined 
by calculating a set of observables {X^l^ } for several spherical nuclei and by 
adjusting the parameters to minimize a generalized defined by |21j : 



where i runs over the set of nuclei, X runs over the set of observables, the 
subscript "exp" indicates the experimental value of the observable, and 
are the relative weights. The weights were chosen to be the expected accu- 
racy for the given observable in a good fit. In practice, a reasonable range of 
weights was tested, and the qualitative conclusions discussed below were al- 
ways reproduced. Some of the considerations relevant to the choice of weights 
are discussed in |S1 15| . 

The relativistic mean-field equations are solved self-consistently for the 
closed-shell nuclei "'^^O, ''"Ca, ''^Ca, ^^Sr, and ^°^Pb, and also in the nuclear 
matter limit. The parameters are then fitted to empirical properties of the 
charge densities, the binding energies, and various splittings between energy 
levels near the Fermi surface using the figure of merit (x^) in H41|l . The full set 
of {X*^*'} comprised 29 calculated and empirical values. When working at the 
highest order of truncation (essentially v = 4), the calculated results are very 
accurate, as we illustrate shortly, but they are too numerous to reproduce 
here [SI E21 ^| ■ Some of the fitted parameters at various levels of truncation 
are given in Tabled 

To simplify the initial discussion, we restrict consideration to infinite nu- 
clear matter. For symmetric matter [N — Z), the energy density through 
order ly — A is given by 



where fcp is the Fermi wavenumber, = 2fcp/37r^, and = ^s^o ~ M — M* 
and W = gvVo are the scaled fields defined earlier in terms of the scalar 
and vector mean fields (j^o and Vq. For readers who are familiar with the 
corresponding result in the Walecka model (Eq. (3.53) in [H]), one sees that 
the MFT nuclear matter energy has been generalized to include additional 
nonlinearities that are not allowed in the renormalizable Walecka model. The 
fields and W are determined by extremization of E. 

What causes the nuclear matter saturation and the relatively small bind- 
ing energy? Let's expand ^mft/Pb from (|42|l in powers of jH) and suppress 




(41) 
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the nonlinear meson terms for clarity: 
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The lowest-order Lorentz scalar and vector contributions (which are propor- 
tional to /Ob) set the scale for the large mean fields 'P and W. This scale is 
consistent with chiral QCD counting rules |H1EII, but these two terms cancel 
almost exactly in the binding energy, leading to an anomalously small remain- 
der. However, they add constructively in the spin-orbit interaction, leading 
to appropriately large spin-orbit splittings in nuclei ^1 121 ■ 

It is important to notice the different behavior of the vector and scalar 
interaction terms in H43() . Whereas the vector interaction enters at only linear 
order in pe, the scalar interaction enters at all orders; moreover, the leading 
scalar term at every order in looks exactly the same, and they all add 
constructively. These terms are precisely what one gets by shifting the nucleon 



Table 1. Parameter sets from fits to finite nuclei, as described in the text. The 
parameters in the lower portion of the table are fitted to the (free) nucleon charge 
and magnetic form factors, and the proton charge e and the anomalous moments 
Ap and An are given their empirical values jH|. 







Wl 


CI 


Ql 


Q2 


Gl 


G2 


ms/M 


2 


0.60305 


0.53874 


0.53735 


0.54268 


0.53963 


0.55410 


ffa/47r 


2 


0.93797 


0.77756 


0.81024 


0.78661 


0.78532 


0.83522 


5v/47r 


2 


1.13652 


0.98486 


1.02125 


0.97202 


0.96512 


1.01560 


Sp/47r 


2 


0.77787 


0.65053 


0.70261 


0.68096 


0.69844 


0.75467 


r?i 


3 




0.29577 






0.07060 


0.64992 




3 




1.6698 


1.6582 


1.7424 


2.2067 


3.2467 


Vp 


3 










-0.2722 


0.3901 


V2 
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-0.96161 


0.10975 
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-8.4836 
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Co 


4 








-1.7750 


3.5249 


2.6416 


Ql 
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1.8549 


1.7234 


Q2 


5 










1.7880 


-1.5798 


/v/4 


3 










0.1079 


0.1734 


/p/4 


3 


0.9332 


1.1159 


1.0332 


1.0660 


1.0393 


0.9619 




4 


-0.38482 


-0.01915 


-0.10689 


0.01181 


0.02844 


-0.09328 




4 


-0.54618 


-0.07120 


-0.26545 


-0.18470 


-0.24992 


-0.45964 
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mass in the nonrelativistic kinetic energy term Sfcp/IOM from M — > M* w 
M — glpj^/ml. These additional, repulsive, velocity-dependent interactions 
reduce the strength of the lowest-order, attractive scalar contribution and 
are crucial for establishing the location of the equilibrium point of nuclear 
matter. Thus the different behavior of the vector and scalar interactions leads 
to large relativistic interaction effects in the nuclear matter energy density. 
In contrast, the relativistic corrections to the kinetic energy (the nonleading 
terms in the first pair of square brackets) are indeed small; this is not where 
the important "relativity" is. 

The critical question is whether the hierarchal organization of interaction 
terms is actually observed. This is illustrated in Fig. where the nuclear 
matter energy /particle is shown as a function of the power of the mean fields, 
which is called 6 in (There are no gradient contributions in nuclear matter 
and (tt) — 0.) The crosses and error bars are estimates based on NDA and 
naturalness, that is, overall coefficients are of order unity. It is clear that each 
successive term in the hierarchy is reduced by roughly a factor of five, which 
implies a value of A w 600 MeV. Thus for any reasonable desired accuracy, 
the Lagrangian can be truncated at a low value of v. In fact, contributions 
to the energy /particle that are smaller than roughly 1 MeV are below the 
level of resolution and can be eliminated in favor of small adjustments of 
the remaining parameters. Derivative terms and other coupling terms will be 
discussed later. 
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Fig. 1. Nuclear matter energy/particle for two QHD parameter sets, one on the 
left (Gl) and one on the right (G2) of the error bars. The power of fields is b = j + £ 
for a term of the form $-'W^ {£ is even). The boxes denote terms with j = 0, the 
circles denote terms with £ — 0, and absolute values are shown. The crosses with 
error bars are estimates based on 13811 . with 1/2 < C < 2. The arrow indicates the 
total binding energy eo = 16.0 MeV. 
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The quality of the fits to finite nuclei and the appropriate level of trun- 
cation is illustrated in Fig. 13 ^H], where the figure of merit is plotted as a 
function of truncation order and of various combinations of terms retained 
in >Cqhd- The full calculations (•) retain all allowed terms at a given level 
of u, while the other two choices keep only the indicated subset. There is 
clearly a great improvement in the fit (more than a factor of 35) in going 
from V = 2 io V = 4,^ but there is no further improvement in going to = 5, 
using the extra interactions contained in H4()|l . Speaking chronologically, the 
V = 2 results show the level of accuracy obtained more than 20 years ago |41| . 
while the i/ = 4 results were obtained seven years ago Moreover, the "^^ 
only" results at = 4 (■) show the state of the situation in the late 1980s, 
as discussed in ^| C] . Recent work JHI shows that the full complement of 
parameters at order i/ = 4 is underdetermined, and that only six or seven are 
determined by this data set, which explains the success of the earlier MFT's 
with a restricted set of parameters 021 ■ We will review this analysis in the 
next section. 

As a further example of the quality of the fits, and for some additional 
historical perspective. Fig. O shows the percent deviation between the cal- 
culated and empirical binding energies for the five closed-shell nuclei listed 
earlier. The results labeled "Walecka" are from those labeled "QMC" 
are from recent quark-meson models of nuclear structure 0S| 0^1 HOI > and 
those labeled "QHD" follow from the present EFT for various parameter sets 
listed in Tabled It is obvious that the modern EFT approach improves the 
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Fig. 2. values for QHD parameter sets, as a function of the level of truncation. 
The power f = 2 corresponds to set Wl, i/ = 3 is for set CI, the v = 4 square 
is for set Ql, the u = 4 diamond is for set Q2, and the circle is for the full set of 
parameters Gl. The v = 5 results include the terms from £5 in 140 l l. 
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Fig. 3. The deviations (in percent) between the calculated binding energy /nucleon 
and the empirical values for five doubly-magic nuclei. The different calculations and 
parameter sets are discussed in the text. For orientation, the Walecka-model results 
are ttnder bound. 



quality of the fits by roughly two orders of magnitude and estabhshes a new 
standard of accuracy that must be attained for any modern approaches to 
nuclear structure to be considered viable. 

Finally, to illustrate the power of the EFT approach to nuclear many- 
body physics, we show two recent results of Huertas The calculations 
use the MFT of -Cqhd with parameters fitted to closed-shell nuclei along the 
"valley of stability", namely, set G2 in Tabled The resulting self-consistent, 
relativistic Kohn-Sham and meson field equations are solved to calculate 
the properties of Sn isotopes out to doubly-magic values of N and Z far 
from stability. These results are therefore true predictions of the EFT, since 
no adjustments were made to previously determined parameters. Figure 0] 
shows the predicted binding energies of the even-even isotopes from ^^qSh to 
""^I^Sn compared with measured experimental results. The theoretical values 
are accurate to better than 1% throughout. (Similar accuracy is obtained for 
parameter set Gl; see Fig. 4 in WB ) 

Figure |S1 shows the predicted ground-state energies, spins, and parities of 
the neighboring single-particle and single-hole nuclei Sb and ^Igln relative 
to "'^ioSn. The energy differences are just the chemical potentials, which should 
be accurately reproduced, according to the discussion of Kohn-Sham theory 
in Sect.O The agreement provides compelling evidence that QHD is indeed 
an EFT for low-energy QCD that can be used to describe nuclear many- 
body physics. This recent work has been extended to semi-magic nuclei with 
N = 28, 50, 82, 126 and Z = 28, 50, 82 in [S]. 
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How can we understand the excellent accuracy of the preceding MFT re- 
sults? As discussed in Sect. 01 the exact energy functional has kinetic-energy 
and Hartree parts (which are combined in a relativistic formulation) plus an 
exchange-correlation functional, which is generally a nonlocal, nonanalytic 
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functional of the densities that contains aU the other many-body, relativis- 
tic, and short-range effects The basic idea behind the relativistic MFT 
(RMFT) is to approximate the functional using an expansion in classical 
meson fields (or nucleon densities) and their derivatives, based on the ob- 
servation that the ratios of these quantities to the nucleon mass are small, 
at least up to moderate density.* The parameters introduced in the expan- 
sion are fitted to experiment, and if we have a systematic way to truncate 
the expansion, the framework is predictive. Moreover, if the RMFT energy 
functional is sufficiently general, it will automatically incorporate effects be- 
yond the Hartree approximation, such as those due to short-range physics 
and many-body correlations. 

But why should we expect an approximate, mean-field energy functional 
to work so well? We observe that while the mean scalar and vector potentials 
# and W are small compared to the nucleon mass, they are large on nuclear 
energy scales [^El- Moreover, as is illustrated in Dirac-Brueckner-Hartree- 
Fock (DBHF) calculations 031123 EH i the scalar and vector potentials (or 
self-energies) are nearly state independent and are almost equal to those 
obtained in the Hartree approximation. Thus the Hartree contributions to 
the energy functional should dominate, and an expansion of the exchange- 
correlation functional in terms of mean fields should be reasonable. This 
"Hartree dominance" also implies that it should be a good approximation 
to associate the single-particle Dirac eigenvalues with the empirical nuclear 
energy levels, at least for states near the Fermi surface p^.^ 

We also observe that the nuclear properties of interest discussed in Sect.Q 
include: 1) nuclear shape properties, such as charge radii and charge densi- 
ties, 2) nuclear binding-energy systematics, and 3) single-particle properties 
such as level spacings and orderings, which reflect spin-orbit splittings and 
shell structure. Since the Kohn-Sham approach is formulated to reproduce 
exactly the ground-state energy and density, and the Hartree contributions 
are expected to dominate the Dirac single-particle potentials, these observ- 
ables are indeed the ones for which meaningful comparisons with experiment 
should be possible. 

As discussed above, an RMFT energy functional of the form in (|42ll . 
extended to include low-order derivatives of the meson fields, successfully 
reproduces these nuclear observables with parameters of natural size (see Ta- 
ble ^1. This justifies a truncation of the energy functional at the first few 
powers of the fields and their derivatives, as is evident from Fig. More- 

^ Since the meson fields are roughly proportional to the nuclear density, and since 
the spatial variations in nuclei are determined by the momentum distributions of 
the valence-nucleon wave functions, this organizational scheme is essentially an 
expansion in kp/M, for fcp corresponding to ordinary nuclear densities. Here the 
nucleon mass M « A is a generic large mass scale characterizing physics beyond 
the Goldstone bosons. 

^ One expects the KS spin-orbit splittings to be more accurate than the absolute 
energy eigenvalues. 
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over, the full complement of parameters is underdetermined, so keeping only 
a subset does not preclude a realistic fit to nuclei. Both the early RMFT cal- 
culations mentioned above and the newer calculations based on chiral EFT 
should be interpreted within the context of this Kohn-Sham approach to 
DFT. 



6 Analysis of Mean-Field-Theory Parameters 

Although we could use meson-nucleon EFT's as done historically, the anal- 
ysis is more transparent with point-coupling theories, which contain only 
nucleon fields in a local Lagrangian. Because of the freedom to perform field 
redefinitions, a general point-coupling Lagrangian is equivalent to a general 
meson-nucleon Lagrangian [HI 122 12I| ■ 

An energy functional of nucleon densities can be constructed by start- 
ing with a general point-coupling effective Lagrangian, consistent with the 
symmetries of QCD, and by evaluating the corresponding one-loop energy 
functional. As discussed above, this approach approximates a general DFT 
functional that incorporates many-body effects beyond the Hartree level when 
the parameters are determined from finite-density data. 

To arrive at a suitable truncation scheme, we again rely on NDA and natu- 
ralness. The two relevant mass scales are and A, and for closed-shell nuclei, 
the energy functional is an expansion in powers of the nucleon scalar, vector, 
isovector- vector, tensor, and isovector-tensor densities, which are defined as 
Ps = (iptp), Pb = (V'V), P3 = liip'^Tsi/j), Si = {ipa^'i)), and ssi = ^{'i{^a°''T3ip), 
respectively, where ip is the nucleon field. 

We can then define scaled densities and their derivatives as [see (|38|) ] 

NDA also provides numerical estimates for the scaled densities that will allow 
us to estimate terms in the energy functional. For example, each additional 
power of Ps is accompanied by a factor of /^A. The ratios of scalar and vector 
densities to this factor at nuclear matter equilibrium density are between 1 /4 
and 1/7 jSTj, which serves as an expansion parameter. Similarly, one can 
anticipate good convergence for gradients of the densities, since the relevant 
scale for derivatives in finite nuclei is the nuclear surface thickness cr, and 
so the dimensionless expansion parameter is 1/Aa < 1/5. The expansion is 
useful because the coefficients have been shown empirically to be natural, 
that is, of order unity pil22|. 

The energy functional is dominated by the isoscalar terms, but we also 
include the isovector and tensor terms for completeness. To nominal order 
ly = A (all densities and gradients count as = 1, except the three- vector 
tensor densities s, which have w 2; see also footnote on p.Ejl, we obtain 
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occ „ 
a 

+ C2Pb - Cd(VpB)' + CiPb - 5iPs(Vp,)2 - 52Ps(Vpb)2 

+ 6P3^ - 6i(Vp3)^ + ?ypP3^Ps + /v Vpb • S + JpVpz ■ St 

+ electromagnetic and higher-order terms , (45) 

where the notation of is used. The sum runs over occupied nucleon states. 

In Fig. 1^ the small symbols show the NDA estimates (with associated 
error bars) for the various energy contributions in ^^O and ^°^Pb. The mag- 
nitudes of energy contributions in l|45(l from two representative RMF point- 
coupling models (i.e., two different parameter sets) are shown as larger un- 
filled symbols (one model on each side of the error bars). These models pro- 
vide very accurate predictions of bulk nuclear properties |22| . The energy 
contributions are determined for each nucleus by making multiple runs while 
varying each parameter slightly around its optimized value, which enables 
us to deduce the logarithmic derivative with respect to each parameter. The 
filled symbols denote the sum of the values for each power of the density. The 
binding energy/nucleon in equilibrium nuclear matter is denoted by eo. 

The two representative point-coupling models validate the isoscalar esti- 
mates (small open squares) , and the resulting hierarchy of isoscalar contribu- 
tions is quite clear. How far down in the hierarchy can we reliably determine 
contributions and their associated parameters? In Fig. 12 the impact of dif- 
ferent truncations of RMF meson-nucleon models is shown by plotting the 
figure of merit against the maximum power of fields. We have also performed 
this test with point-coupling models, with similar conclusions. The "full" 
models (which include all nonredundant terms at a given order) show that 
one needs to go to the fourth power of the fields to get the best fits, but going 
further yields no improvement. Analogous behavior is found for RMF point- 
coupling models with powers of densities replacing powers of fields i22j . Thus 
contributions to the energy/particle at the level of roughly 1 MeV are at the 
limit of resolution. Fifth-order isoscalar contributions to the energy /particle, 
which are predicted to be less than 1 MeV, are simply not determined by the 
optimization. 

The variation in coefficient values in H45|) provides a measure of how well 
the parameters are actually determined by the data. Figure [3 shows the 
seven coefficients of isoscalar non-gradient terms from four point-coupling 
models. Note that all coefficients are natural, i.e., of order unity. However, 
the spread in coefficient values is significant and does not correspond well to 
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the power-counting order. We conclude that different linear combinations of 
the coefficients must be considered to draw reliable conclusions about how 
many are determined by the data. 
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Fig. 6. Contributions to the energy /particle in ^ O and Pb determined by log- 
arithmic derivatives with respect to the model parameters (see text) for two RMF 
point-coupling models |22| . Absolute values are shown. The filled symbols are net 
values. The small symbols indicate estimates based on NDA |10) . with the error 
bars corresponding to natural coefficients from 1/2 to 2. 



Can we find a more systematic power counting scheme? The similar size 
of the scalar density ps and the vector density ps suggests that we count 
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Fig. 7. Coefficients for four accurately fit RMF point-coupling models from |22) . 
Each model is represented by a different shape, and the shading shows the type of 
term (scalar, vector, or mixed). 



instead powers of p+ = (ps + Pb)/2 and p_ = (ps — Pb)/2. The corresponding 
"improved" coefRcients are listed in Table |21 [see (|45f) ]. The spread in these 
coefficients for four RMF point-coupling models from 22 are shown in Fig.|Sl 
The terms are organized according to the powers of p+ and with p- 
scahng as p^J_^ . 

The leading orders are very well determined, with a systematic^ increase 
in uncertainty. Even the sign is undetermined for the parameter 124, which 
is shown with unfilled symbols, but the next parameter (i^e) appears to 
be reasonably well determined. Higher-order terms are not determined by 
the optimizations. Deduced values and uncertainties based on this sample of 
models are given in Tabled We see that of the seven isoscalar non-gradient 
parameters in H45|) . four linear combinations are clearly determined by bulk 
nuclear observables, with probably a fifth combination as well. 



Table 2. Improved coefficients for point-coupling RMF models (Table VI of |22|). 
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Fig. 8. Improved coefficients for tlie same four models as in Fig. |7| The "order" is 
determined by counting powers of p+ and p_ . 



The isovector terms appear only on the graphs for ^°®Pb. The factor 
{N — Z)/2B, which is only 10% even for Pb, severely limits the sensitivity to 
isovector terms (especially since the factor must appear with even powers). 
The magnitude of the leading isovector term (oc £,2pi ) is comparable to the 
fourth-order isoscalar term, which is at the limit of what can be determined 
reliably from fitting the binding energy. We conclude that only one isovector 
parameter is determined by the bulk observables. See the contribution to 
these lectures by R. Furnstahl for further discussion of the isovector terms. 

The energy estimates for isoscalar, tensor terms imply that only one pa- 
rameter, at best, can be determined. Higher-spin terms, which will require 
more gradients and have smaller average densities, are not at all constrained. 
The tensor terms are interesting because a fraction of the spin-orbit force can 
be generated by including an isoscalar, tensor coupling of the vector field to 
the nucleon |14) . Nevertheless, the spin-orbit potential arises predominantly 
from the large scalar and vector fields; attributing more than one-third of the 
potential to the tensor coupling produces unrealistic surface systematics |58| . 

Finally, the gradient terms follow the same pattern in the energy: the 
leading term is barely above the limit of resolution. In fact, there are two 
isoscalar gradient terms at leading order (scalar and vector), but only their 
sum is well determined. The sum of the subleading-order contributions almost 
vanishes. Thus we conclude that only one gradient parameter is determined. 

The handful of parameters that are well determined by the usual bulk 
nuclear observables (binding energies, charge density distributions, and spin- 
orbit splittings in doubly magic nuclei) can be associated with an equal num- 
ber of nuclear properties and general features of RMFT's. In particular, 

1.) Two isoscalar non-gradient parameters are very well determined. These 
correspond to the highly constrained values for the equilibrium density 
(fcp w 1.30 ± 0.01 fm"^) and binding energy (16.0 ± 0.1 MeV) of nuclear 
matter. 
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2. ) An additional isoscalar constraint is that AI* « 0.61 ± 0.03, if the 

isoscalar, tensor term is set to zero. This range ensures an accurate 
reproduction of spin-orbit spUttings in finite nuclei. Small increases in 
M* without changing the splittings can be accomplished by including an 
isoscalar, tensor term; an analysis using a simple local-density approxi- 
mation is discussed in |14| . 

3. ) A fourth isoscalar constraint comes from the nuclear matter compress- 

ibility. The constraint is much weaker, in the range of X « 250 ±50 MeV. 

4. ) The possibility of a fifth isoscalar constraint has been considered by 

Gmuca jSHl, who argued that separate scalar and vector fourth-order 
terms were necessary to tune the density dependence of the scalar and 
vector parts of the baryon self-energy. This would correspond to con- 
straining [2q from Table |21 Moreover, some form of isoscalar nonlinear 
vector interaction is needed to soften the high-density equation of state 
to be consistent with observed neutron star masses j53) . 

5. ) Since only one isoscalar gradient parameter is determined, it is not useful 

to allow the scalar and vector masses (or their equivalents in a point- 
coupling theory) to vary independently. Thus it is convenient to fix the 
vector mass at a natural size, such as the experimental mass for the lv. 
A scalar mass of 500 ± 20 MeV is then required. 

6. ) The one isovector parameter can be fixed by the surface-corrected volume 

symmetry energy jgO, , which falls in the range 34 ± 4 MeV 0^ . Since no 
isovector gradient is determined, setting the isovector, vector meson mass 
to the experimental p meson mass is adequate. 

To ensure a reasonable (if not optimal) description of finite nuclei, it is 
sufficient to reproduce the nuclear matter properties given above, but note 
that all properties must be satisfied, and the resulting parameters must be 
natural. One cannot justify the underlying physics of a model if it repro- 
duces only a subset of the nuclear calibration data. 

7 Weak Nuclear Currents 

A desirable theory of nuclear currents should satisfy the following three cri- 
teria: 

• It should use the same degrees of freedom to describe the currents and 
the strong-interaction dynamics. 

• It should satisfy the same internal symmetries, both discrete and contin- 
uous, as the underlying theory of QCD. 

A further caution is that there are many correlations among these properties, so 
that the allowed ranges should not be considered to be independent. 
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• Its parameters can be calibrated using strong-interaction phenomena, like 
ttN scattering and the properties of finite nuclei. This is especially impor- 
tant in eft's, as they contain all (non-redundant) interaction terms that 
are consistent with the underlying symmetries (611 [5] . 

The QHD framework described so far embodies these three desirable features. 

The weak currents arise from Noether's theorem applied directly to £qhd 
and contain the pion field to all orders. The leading-order (in i/) vector and 
axial-vector currents are given by (a is the isospin index) 

f2 1 _ 

+ I QaN^^I, [^T'^e - eV^e] N , (46) 

f2 1 

= -i ^ Tr {r° (C/9^[/^ - U^d^'U) ] - - N-f^ [^r^^t - £}t°-£] N 

- \ QaN^^I, [ir'^e + eV"^] N . (47) 

As shown in |S2| , in the presence of an external axial- vector source, the scat- 
tering amplitudes constructed from these currents satisfy CVC, PCAC (when 
rriTr ^ 0) and the Goldberger-Treiman relation (with t;^ ^ 1) automatically. 
Moreover, the chiral charges Q"" and Qg derived from these currents satisfy 
the familiar chiral charge algebra to all orders in the pion field. In 62 , these 
currents are used to study beta decay in ^■^^'^^■^Sn. 



8 Summary 

In this talk I discussed recent progress in Lorentz-covariant quantum field 
theories of the nuclear many-body problem, often called quantum hadrody- 
namics (QHD). QHD is a local, nonrenormalizable, effective Lagrangian field 
theory with baryons and mesons as the generalized coordinates (fields). An 
effective Lagrangian consists of known long-range interactions constrained by 
symmetries and a complete, non-redundant set of short-range interactions. 
By simply looking at the spectra of massive nuclei, it is obvious that some 
relativistic effects must be important in nuclei; thus, it is most convenient to 
use a Lorentz-covariant theory. 

The effective field theory studied here contains nucleons, pions, isoscalar 
scalar [a) and vector (cj) fields, and isovector vector {p) fields. The heavy 
mesons are introduced as collective, effective degrees of freedom to simplify 
the description of the medium- and short-range nucleon-nucleon interaction 
and to conveniently parametrize ground-state expectation values of nucleon 
bilinears, which are important for the description of bulk nuclear properties. 
The QHD theory exhibits a nonlinear realization of spontaneously broken 
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SU(2)^ X SU(2)^ chiral symmetry and has three desirable features: it uses 
the same degrees of freedom to describe the nuclear currents and the strong- 
interaction dynamics, it satisfies the symmetries of the underlying theory of 
QCD, and its parameters can be calibrated using strong- interaction phenom- 
ena. Moreover, the electromagnetic structure of the nucleon can be included 
straightforwardly in a derivative expansion of the fields. 

Although the QHD Lagrangian in principle contains an infinite number 
of terms, naive dimensional analysis and naturalness allow one to identify 
suitable expansion parameters and to estimate the sizes of various terms in 
the Lagrangian. Thus, for any desired accuracy, the Lagrangian can be trun- 
cated to a finite number of terms. In particular, for normal nuclear systems, 
it is possible to expand the QHD effective Lagrangian systematically in pow- 
ers of the meson fields (and their derivatives) and to truncate the expansion 
reliably after the first few orders. 

Using density functional theory, I showed that the mean-field approxima- 
tion produces an energy functional whose parameters can be determined by 
fitting bulk and single-particle properties of nuclei. The framework of Kohn 
Sham theory allows the ground state to be constructed from (quasi)particle 
orbitals with unit occupation number. Since the mean-field energy functional 
is a good approximation to the exact energy functional over the relevant range 
of density, it is possible to reproduce nuclear densities, binding energies, and 
single-particle spectra near the Fermi surface very accurately. Because the 
parameters are fitted to nuclear properties, the energy functional implic- 
itly contains effects that go beyond a simple Hartree approximation, such as 
short-range physics, hadron substructure, and many-body correlations. 

The numerical parameters of QHD were studied using an effective, point- 
coupling Lagrangian that contains nucleon fields only. Because of the free- 
dom to redefine the fields (or coordinates), a general point-coupling theory 
is equivalent to a general baryon meson theory. By examining the contribu- 
tions to the energy/nucleon in doubly magic nuclei, it was found that only a 
small number of parameters (roughly seven) can be calibrated by the nuclear 
data input. New ways to calibrate additional parameters will play an impor- 
tant role in the construction of the next generation of QHD Lagrangians, as 
discussed by R. Furnstahl elsewhere in this volume. Finally, the weak vector 
and axial- vector currents in the QHD framework were discussed. 

Nuclear physics is the study of strongly interacting hadronic matter, and 
the only consistent theoretical framework we have for describing such a rel- 
ativistic. interacting, quantum-mechanical, many-body system is relativistic 
quantum field theory based on a local Lagrangian density. Although QCD of 
quarks and gluons provides the basic underlying theory, Lagrangians com- 
prised of hadronic degrees of freedom (QHD) provide the most efficient de- 
scription of the physics in the strong-coupling domain. In the modern effective 
field theory perspective of QCD, one incorporates only the underlying sym- 
metries of QCD into the QHD Lagrangian. By interpreting the mean-field 
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approximation in the context of density functional theory, one can under- 
stand the numerous successes in the QHD description of nuclear properties. 
Nevertheless, finding an efficient, tractable, nonperturbative way to match 
the QCD Lagrangian to the long-range, strong-coupling, effective field the- 
ory of QHD is still a major goal for the future. 

Acknowledgments 

I am grateful to Dick Furnstahl and Dirk Walecka for useful comments on a 
draft of this manuscript. This work was supported in part by the US Depart- 
ment of Energy under Contract No. DE-FG02-87ER40365. 

References 

1. A. L. Fetter and J. D. Walecka: Quantum Theory of Many-Particle Systems 
(McGraw-Hill, New York 1971) chaps 11 and 15. Reprinted by Dover Publi- 
cations, Mineola, New York 2003 

2. J. D. Walecka: Theoretical Nuclear and Subnuclear Physics (Oxford Univ. 
Press, New York 1995) part I 

3. B. D. Serot and J. D. Walecka: Adv. Nucl. Phys. 16, 1 (1986), and references 
therein 

4. J. D. Bjorken and S. D. Drell: Relativistic Quantum Mechanics (McGraw-Hill, 
New York 1964, 1998) Fig. 4- 2 

5. M. G. Mayer and J. H. D. Jensen: Elementary Theory of Nuclear Shell Struc- 
ture (Wiley, New York 1955) p 58 

6. J. D. Walecka: Ann. Phys. (New York) 83, 491 (1974) 

7. B. D. Serot: Rep. Prog. Phys. 55, 1855 (1992) 

8. R. J. Furnstahl, B. D. Serot, and H.-B. Tang: Nucl. Phys. A615, 441 (1997); 
A640, 505 (1998) (E) 

9. B. D. Serot and J. D. Walecka: Int. J. Mod. Phys. E 6, 515 (1997), and refer- 
ences therein 

10. R. J. Furnstahl and B. D. Serot: Nucl. Phys. A671, 447 (2000) 

11. R. J. Furnstahl and B. D. Serot: Nucl. Phys. A673, 298 (2000) 

12. R. J. Furnstahl and B. D. Serot: Comm. Mod. Phys. 2, A23 (2000) 

13. B. D. Serot and J. D. Walecka, in: 150 Years of Quantum Many-Body The- 
ory, ed by R. F. Bishop, K. A. Gernoth, and N. R. Walet (World Scientific, 
Singapore 2001) pp 203-211 

14. R. J. Furnstahl, J. J. Rusnak, and B. D. Serot: Nucl. Phys. A632, 607 (1998) 

15. T. D. Cohen, R. J. Furnstahl, and D. K. Griegel: Phys. Rev. Lett. 67, 961 
(1991) 

16. R. J. Furnstahl, D. K. Griegel, and T. D. Cohen: Phys. Rev. C 46, 1507 (1992) 

17. R. J. Furnstahl, H.-B. Tang, and B. D. Serot: Phys. Rev. C 52, 1368 (1995) 

18. H.-B. Tang: hep-ph/9607436, unpublished 

19. P. J. Ellis and H.-B. Tang: Phys. Rev. C 57, 3356 (1998) 

20. T. Becher and H. Leutwyler: Eur. Phys. J. C9, 643 (1999) 



32 Brian D. Serot 



21. B. A. Nikolaus, T. Hoch, and D. G. Madland: Phys. Rev. C 46, 1757 (1992) 

22. ,]. ,]. Rusnak and R. J. Furnstahl: Nucl. Phys. A627, 495 (1997) 

23. H.-B. Tang and P. J. Ellis: Phys. Lett. B387, 9 (1996) 

24. H. Gcorgi and A. Manohar: Nucl. Phys. B234, 189 (1984) 

25. H. Georgi: Phys. Lett. B298, 187 (1993) 

26. R. J. Furnstahl, B. D. Serot, and H.-B. Tang: Nucl. Phys. A598, 539 (1996) 

27. R. J. Furnstahl, H.-W. Hammer, and N. Tirfcssa: Nucl. Phys. A689, 846 (2001) 

28. B. D. Serot: Phys. Lett. 86B, 146 (1979); 87B, 403 (1979) (E) 

29. S. Coleman, J. Wess, and B. Zumino: Phys. Rev. 177, 2239 (1969) 

30. C. G. Callan, Jr., S. Coleman, J. Wess, and B. Zumino: Phys. Rev. 177, 2247 
(1969) 

31. U. van Kolck: Prog. Part. Nucl. Phys. 43, 337 (1999) 

32. S. M. Ananyan, B. D. Serot, and J. D. Walecka: Phys. Rev. C 66, 055502 
(2002) 

33. R. M. Dreizler and E. K. U. Gross: Density Functional Theory (Springer, Berlin 
Heidelberg New York 1990) 

34. L. H. Thomas: Proc. Cambridge Phil. Soc. 23, 542 (1927) 

35. E. Fermi: Z. Phys. 48, 73 (1928) 

36. N. Argaman and G. Makov: Am. J. Phys. 68, 69 (2000) 

37. W. Kohn: Rev. Mod. Phys. 71, 1253 (1999) 

38. W. Kohn and L. J. Sham: Phys. Rev. A140, 1133 (1965) 

39. A. H. MacDonald and S. H. Vosko: J. Phys. C 12, 2977 (1979) 

40. C. O. Almbladh and U. von Barth: Phys. Rev. B 31, 3231 (1985) 

41. C. J. Horowitz and B. D. Serot: Nucl. Phys. A368, 503 (1981) 

42. P. Ring: Prog. Part. Nucl. Phys. 37, 193 (1996), and references therein 

43. C. J. Horowitz and B. D. Serot: Nucl. Phys. A464, 613 (1987); A473, 760 
(1987) (E) 

44. J. L. Friar, D. G. Madland, and B. W. Lynn: Phys. Rev. C 53, 3085 (1996) 

45. S. Weinberg: Phys. Lett. B251, 288 (1990) 

46. W. H. Furry: Phys. Rev. 50, 784 (1936) 

47. R. J. Furnstahl, C. E. Price, and G. E. Walker: Phys. Rev. C 36, 2590 (1987) 

48. K. Saito and A. W. Thomas: Phys. Lett. B327, 9 (1994) 

49. K. Saito, K. Tsushima, and A. W. Thomas: Nucl. Phys. A609, 339 (1996) 

50. K. Saito, K. Tsushima, and A. W. Thomas: Phys. Rev. C 55, 2637 (1997) 

51. M. A. Hucrtas: Phys. Rev. C 66, 024318 (2002); 67, 019901 (2003) (E) 

52. M. A. Huertas: nucl-th/0304016. 

53. H. Miiller and B. D. Serot: Nucl. Phys. A606, 508 (1996) 

54. A. R. Bodrncr: Nucl. Phys. A526, 703 (1991) 

55. B. ter Haar and R. Malfiict: Phys. Rep. 149, 207 (1987) 

56. R. Machleidt: Adv. Nucl. Phys. 19, 189 (1989) 

57. J. L. Friar: Few Body Sys. 99, 1 (1996) 

58. G. Hua, T. v. Chossy, and W. Stocker: Phys. Rev. C 61, 014307 (2000) 

59. S. Gmuca: Z. Phys. A342, 387 (1992); Nucl. Phys. A547, 447 (1992) 

60. P. A. Seeger and W. M. Howard: Nucl. Phys. A238, 491 (1975) 

61. S. Weinberg: The Quantum Theory of Fields, vol. I: Foundations (Cambridge 
Univ. Press, Cambridge, UK 1995) 

62. M. A. Huertas: Acta Phys. Polon. B 34, 4269 (2003) 



